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The pyrochlore spin ice compounds Dy2Ti207 and Ho2Ti207 are well described by classical Ising models 
down to very low temperatures. Given the empirical success of this description, the question of the importance 
of quantum effects in these materials has been mostly ignored. We show that the common wisdom that the 
strictly Ising moments of Dy^^ and Ho^+ imply strictly Ising interactions is too naive; a more complex argument 
is needed to explain the close agreement between theory and experiment. From a microscopic picture of the 
interactions in rare-earth oxides, we show that the high-rank multipolar interactions needed to induce quantum 
effects are generated only very weakly by super-exchange. Using this framework, we formulate an estimate of 
the scale of quantum effects in Ho2Ti207 and Dy2Ti207, finding it to be well below experimentally relevant 
temperatures. We discuss the implications of these results for realizing quantum spin ice in other materials. 


Introduction: Spin ice has proven to be one of the more 
fruitful marriages of theoretical and experimental condensed 
matter physics [1-5]. This cooperative paramagnetic [6] 
phase is a magnetic analogue of common water ice [ 1 , 2] , with 
proton displacements mapped to magnetic moments pointing 
in or out of the comer-shared tetrahedra of the pyrochlore 
lattice [7]. Spin ice displays an exponential number of low- 
energy states, and thus an associated extensive residual en¬ 
tropy [8, 9]. This manifold is characterized by the ice rules 
specifying that on each tetrahedron two spins must point in, 
and two must point out. Genetically referred to as a “Coulomb 
phase” [10, 11], the spin ice state harbors a rich phenomenol¬ 
ogy such as bow-tie shaped singularities (pinch points) in 
the magnetic structure factor [11, 12] and gapped low-energy 
excitations which provide a condensed matter realization of 
“magnetic monopoles” [4, 10]. 

The two textbook materials that realize this spin ice physics 
are the pyrochlore rare-earth titanates [3] Ho2Ti207 [7] and 
Dy2Ti207 [8]. Significant evidence has accumulated that 
their magnetic and thermodynamic properties [10, 13-15] 
can be quantitatively described by a classical model that in¬ 
cludes nearest-neighbor exchange and long-range magneto¬ 
static dipole-dipole interactions, both of purely Ising type. 
Amendments incorporating Ising interactions beyond nearest- 
neighbors have been proposed to account for various fine 
details in the thermodynamic and magnetic behavior of 
Dy2Ti207 [14, 16, 17]. The combination of Ising exchanges 
and long-range dipolar interactions is expected to lift the de¬ 
generacy of the ice manifold and release the residual entropy 
at very low temperatures [18-21]. Experimentally, such a 
transition has yet to be observed [22], with further significant 
deviations [23] from predictions raising questions to the com¬ 
pleteness of the classical dipolar Ising spin-ice model [16]. In 
particular, one could speculate that the observed rise of the 
magnetic specific heat below ~ 500 mK [23] is an indication 
that quantum effects are becoming important [24, 25]. 

This classical Ising description [13, 15, 26] of the interac¬ 
tions has framed the theoretical and experimental perspectives 
[1-5] on these materials dating back to their initial discov¬ 
ery [7, 8]. Beyond its empirical successes, this mind-set [13- 


21, 26] is rooted in the single-ion magnetic properties of Dy^^ 
[22, 27-29] and Ho^^ [26, 30-32], specifically in the Ising na¬ 
ture of the crystal electric field (CEF) ground doublet. A con¬ 
sequence of this strict Ising anisotropy is that the transverse 
components of the total angular momentum, J*, vanish be¬ 
tween states of the CEF doublet. This is sufficient to explain 
the Ising form of any bilinear, anisotropic exchange interac¬ 
tions ~ Vj, including the long-range dipolar interactions 
[2]. The theoretical basis of the classical Ising description of 
spin ices then appears to rest on the implicit assumption that 
interactions between the rare-earth ions takes such a bilinear 
form. This common assumption is in fact incorrect, due the 
large spin-orbit coupling in rare-earth ions, strong multipolar 
interactions between the momenta are generated when the 
microscopic super-exchange is downfolded into the free-ion 
manifold [33]. Indeed, this argument has been invoked 
[34] to show that quadrupolar interactions may not be neg¬ 
ligible in Pr-based pyrochlore oxides with strictly Ising mo¬ 
ments, such as Pr2Sn207 [35] and Pr2Zr207 [36]. Further 
compounding the problem is the possibility that even moder¬ 
ate transverse couplings could preempt the ordering expected 
from the dipolar interactions. The existence of other spin ice 
compounds, such as Ho2Sn207 [28, 37], Dy2Sn207 [28, 38], 
Ho2Ge207 [32, 39] and Dy2Sn207 [29, 39] makes an explana¬ 
tion through some accidental fine-tuning unlikely. An unan¬ 
swered puzzle thus lies in the identification of Dy2Ti207 and 
Ho2Ti207 as spin ices: Why are these compounds so well de¬ 
scribed by a classical Ising model down to low temperatures? 
What principles, if any, constrain the scale of the multipolar 
interactions and suppress significant quantum fluctuations? 

In this letter, we explain why quantum effects in Dy- and 
Ho-based spin ice materials are small, providing approxi¬ 
mate upper bounds on their size. First, we outline possible 
sources of exchange interactions and show that only magnetic 
dipole-dipole (MDD) and super-exchange between the rare- 
earth ions are potentially non-negligible at experimentally rel¬ 
evant temperatures. As the MDD interaction is purely Ising 
when projected into the ground CEF doublet, we focus on 
super-exchange. Using a microscopic model of the oxygen 
charge transfer processes, we find that the resulting multipolar 
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inter-ionic couplings arising from super-exchange are strongly 
suppressed beyond rank seven [40]. Because of the composi¬ 
tion of the ground doublet, operators of rank seven or less 
can only connect its sub-leading spectral components and are 
thus highly suppressed. From experimental constraints on the 
size of these sub-leading components, we formulate estimates 
of the magnitude of the transverse (non-Ising) couplings in 
the low-energy theory, estimating them to be roughly two 
orders of magnitude smaller than the nearest-neighbor Ising 
coupling. To make this more concrete, we carry out a model 
calculation of the super-exchange interaction in Dy 2 Ti 207 and 
Ho 2 Ti 207 . This calculation corroborates our argument and 
provides explicit estimates for the size of quantum effects in 
the Ho 2 Ti 207 and Dy 2 Ti 207 spin ices. 

Since the arguments we follow for Dy 2 Ti 207 and Ho 2 Ti 207 
are similar, we focus our presentation on Dy 2 Ti 207 and leave 
a detailed exposition of the results for Ho 2 Ti 207 to the Sup¬ 
plemental Material [41]. 

Pseudo-spin model: Our ultimate goal is to understand the 
interactions acting between the low energy states of the CEF 
manifold. In Dy 2 Ti 207 this is a Kramers doublet built from 
states in the ®Hi 5/2 manifold of the 4/^ configuration of Dy^^. 
For a realistic CEF potential [30, 42], the ground state is a 
doublet with a gap of ~ 30 meV to the first excited level. 
This doublet transforms in the F 5 © Fg representation of the 
site symmetry group and has Ising-like character [43], as 
shown below. One writes the doublet states as 

1+) = Q'|±y)±di I±|)-d2|±|)+d3| + i)+^4| + |)±55 | + t)> 

( 1 ) 

where \M) = 115/2, M) are are eigenstates of the and F op¬ 
erators of the ®Hi 5/2 manifold. We choose the doublet basis so 
that ^ 2^3 -l>6\6n + Sadg = 0 to make the dipole operator di¬ 
agonal. The matrix elements of the magnetic dipole operators 
then take the form 

<(r| = 0 , <(r| !>') = Acrd^^,, ( 2 ) 

where cr - ± and the d-factor is 

d = y - 3 ((5^ -H 2d^ H- 2>5l + A6l + 55^) . (3) 

Including the Dy^^ Lande factor g - 4/3, the magnetic mo¬ 
ment yu = gpB.A is found experimentally to be very close to 
maximal. Estimates from the saturation magnetization and 
high-temperature susceptibility [22, 27] give p = (10+0.1) ps- 
This constrains the |±) states in Eq. (1) to be predominantly 
1+15/2) [ 8 , 30]. We thus have an upper bound of 6^ < 0.025/n 
given + 0.1 ps error bars. 

Since the CEE energy scale is large, it is sufficient to first 
consider only a bare projection of the microscopic Dy^’^-Dy^^ 
interactions acting within the full ®Hi 5/2 manifold into the 
ground doublet. Such a model can be expressed in terms of 
the components of a pseudo-spin 1 /2 operator 

5* = |+)<+|, i(|+)<+|-|-)<-|). (4) 


The symmetry of the pyrochlore lattice constrains the form 
of any such projected Hamiltonian. The nearest-neighbor ef¬ 
fective exchange model for Dy 2 Ti 207 is restricted to the form 
[43] 


(ij) 






( 5 ) 

Since the are not bond dependent, the coupling can 
be removed by a global rotation of the pseudo-spin operators 
about the y axis [43], yielding 


(ij) 


We will find below that the net J^z, including both the con¬ 
tribution from microscopic inter-ionic super-exchange and the 
nearest-neighbor MDD contribution, is the dominant coupling 
and that it is positive, as originally found empirically in Ref. 
[13]. At leading order, the rotated couplings are 

j ~ j _ — J - J J ~ J + — (li 

'fX 'fXX J •) y yy^ z zz ^ j • W/ 

'^ZZ '^ZZ 

This rotation affects the relationship between the components 
of the dipole moment and those of the pseudo-spin 1/2 at 
OiJxzlJzz)- We characterize deviations from the Ising limit us¬ 
ing the nearest-neighbor transverse scale J± = {Jx+Jy)IA [43]. 
Our goal is to estimate the size of J± and thus the strength 
of the quantum effects within the low-energy pseudo-spin 1/2 
model. 

Interactions in rare-earth pyrochlores: To understand how 
the microscopic inter-ionic coupling mechanisms generate 
we need to consider higher energy physics and a descrip¬ 
tion beyond the lowest CEE doublet. At the atomic level, 
there is a hierarchy of energy scales for the rare-earth ions: 
the Coulomb interaction being largest, followed by spin-orbit 
coupling and the CEE potential. The Coulomb interaction 
plays two roles: it provides the repulsion that separates the 
different 4/" charge states and splits each 4/" manifold in ap¬ 
proximate fixed orbital and spin angular momenta [44]. 
Spin-orbit coupling then lifts the degeneracy of the states with 
the same total angular momentum J, giving the terms 
Einally, the single-ion spherical symmetry of the ion is broken 
by the CEE and the manifolds are further split into CEE 
multiple ts. 

We thus move one abstraction level up from the pure ground 
doublet, considering interactions within the ®Hi 5/2 manifold. 
Due to the spherical symmetry of the free ions, it is useful to 
categorize interactions and operators in terms of their rank as 
spherical tensors. Due to the large value J = 15/2 in the ®Hi 5/2 
manifold of Dy 2 Ti 207 , multipolar interactions up to rank-15 
are possible [33]. Interactions at this level will typically be 
on the order of ~ IK, as inferred from the Curie-Weiss scale 
[7, 8 ]. Since the CEE energy scale is of order 300 K [30], 
it is sufficient to consider only the projection into the ground 
doublet. Perturbative corrections are then expected at second 
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order [45] with scale of ~ (1 K)V300 K ~ 3 mK - 5 mK and 
thus can be neglected. There are several sources of interac¬ 
tions between rare-earth ions: electro- and magnetostatic in¬ 
teractions, super-exchange, direct exchange, lattice-mediated 
exchange, etc. Of these the MDD and super-exchange inter¬ 
actions are found to be the most significant [41]. 

The long-range MDD interaction features prominently in 
classical spin ices [13-15, 26]. When projected into the 
ground doublet, it takes the form 


t<j 


/ r 


/ ^nn 1 r„ „ /„ „ \ 

C ' M 



\nj) L V 2 



where is the nearest-neighbor distance between the rare- 
earth atoms, rij = n - rj is the displacement vector between 
the two sites and zi is the quantization axis along the lo¬ 
cal cubic [111] direction at site f]. The strength of the cou¬ 
pling is characterized by The largest 

piece is the nearest-neighbor contribution [16], with coeffi¬ 
cient D s 5^/3 which is ~ 8.9 KinDy 2 Ti 207 [46]. To obtain 
a quantitative agreement with experiments, corrections to the 
Ising couplings must be included. The nearest-neighbour cor¬ 
rection J is the largest, estimated to be ~ 4.96 K in Dy 2 Ti 207 
[13], later rehned to ~ 4.55 K [16]. This coupling has usu¬ 
ally been attributed to super-exchange processes between the 
rare-earth ions and provides an indication of their overall scale 
[47]. We are thus confronted with the question: with such a 
large scale, why does super-exchange not generate significant 
transverse couplings J±1 

Oxygen-mediated super-exchange: To answer this ques¬ 
tion, we consider super-exchange interactions proceeding 
through the oxygen atoms that surround each rare-earth ion. 
To simplify our notation, we single out one exchange path 
consisting of two rare-earths and one oxygen. For the rare- 
earth sites, we define the creation operators and while 
for the intermediate oxygen site we use A combined spin- 
orbital index a = {m, cr) where m and cr label the orbital and 
spin quantum numbers is used throughout. Due to the local¬ 
ized nature of the rare-earth ions, we start with the atomic 
Hamiltonian at each site 


Ho = Hf,i + Hf2 + Hp, (9) 

where Hj i and Hf 2 are the atomic Hamiltonian for the two 
rare-earth ions and Hp for the oxygen site. On the oxygen 
site, we consider only the atomic potential Apj and the cost to 
place two holes together on the same oxygen. Up. We will not 
invoke the details of the rare-earth atomic Hamiltonian, aside 
from selecting the ®Hi 5/2 manifold of the 4/® configuration. 
We perturb Hq with the hybridization terms 

+ l2,al3flPfi + h.C.] , (10) 

a/3 

that represent electron hopping between the orbitals of the 
rare-earth and oxygen ions. Super-exchange interactions are 
generated at fourth-order in perturbation theory in V [34]. For 


the details of the charge-transfer processes, see the Supple¬ 
mental Material [41]. 

To simplify the resolvents that appear in the perturbative 
expansion [41, 48], we follow Ref. [34] and make the com¬ 
mon approximation that only the charging energies - 

are signihcant. This neglects all other intra- 
atomic splittings on the rare-earth sites, such as those due to 
Hund’s coupling or spin-orbit coupling. The energies U/ are 
expected to be of order 5 - 10 eV in 4/ compounds [49]. 
Within this approximation, the effective Hamiltonian can be 
written [34] 

Heff = X (11) 

Q-yS/iV 


where P, projects into the ground state manifold of Hf j at site 
i. The interaction matrix I is defined as [41] 


r 


a(3fiv _ 
12 


(U/ + Apf)2 [ 2U} + Up + 2Apf 


Ud + 2 Au 
\ 




2Uj + Up-\- 2Apf 


Ml [tA\ 


( 12 ) 


Here we dropped constants and single site terms that serve 
only to renormalize the on-site single-ion Hamiltonian. 

Generically, the interactions I are complicated due to the 
states that make up the ®Hi 5/2 ground state manifold of 
Dy2+. However, a strong constraint arises from the one- 
electron form flf/j taken by the operators at each rare-earth 
site. When projected into the ®Hi 5/2 manifold, these operators 
can generate only a small subset of all possible interactions 
between the J = 15/2 degrees of freedom of Dy^^. This can 
be seen explicitly by considering how these operators trans¬ 
form under rotations. Since each / electron carries a total 
angular momentum of J = 5/2 or J = 7/2, the operators f/fp 
are multipoles with rank ranging from rank-0 (5/2 - 5/2 or 
7/2 - 7/2) up to rank-7 (7/2 -H 7/2). Hence, we find that only 
multipolar interactions up to rank-7 are generated by super¬ 
exchange. Because of the spherical symmetry of the free ions; 
the complicated intra-atomic interactions are rank-0 operators 
and thus do not change the rank counting above [41]. 

Including the CEF splittings in the resolvents only changes 
this result slightly. Given the small overall splitting A < 
100 meV induced by the CEF [30, 42], we can treat this as 
a perturbation to Hq along with V. Since the CEF potential is 
a one-electron operator and is time-reversal invariant, it con¬ 
tains only operators up to rank-6 for the same reasons as dis¬ 
cussed above [50]. Each inclusion of the CEE operator can 
thus increase the rank by 6, with up to rank-13 operators at 
order ~ operators up to the maximal rank-15 at 

order ~ (A/t/^)^. The rank-15 operators that directly link the 
leading 1 + 15/2) components are thus strongly suppressed by a 
factor of (A/t/p2 10^“*. We can thus ignore these small cor¬ 

rections and can consider only interactions of rank-7 or lower. 

General argument: Returning to a more general physi¬ 
cal picture, the above results allow us to explain why quan¬ 
tum effects are small in the classical spin ices Dy 2 Ti 207 and 
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Ho 2 Ti 207 . Given the dominant 1+15/2) or |+8) CEF ground 
states in and Ho^^ [30], a large transverse scale J± must 
originate from very high-rank multipole interactions. These 
multipole interactions are only generated weakly from super¬ 
exchange contributions, as they are much higher than rank-7. 
Due to this suppression, the leading contributions to the trans¬ 
verse couplings must come from the sub-dominant spectral 
components of the CEF ground doublet, the components 
given in Eq. (1) for Dy 2 Ti 207 . When these are included, gen¬ 
eration of by super-exchange becomes feasible. As a rough 
estimate, if J is the typical scale of contributions from super¬ 
exchange then we expect 

J± ~ (^) I Jz ~ -J, (13) 


where d represents the maximum size of the sub-dominant 
components of the CEF doublet. Each transverse pseudo-spin 
operator 5* contributes one factor of (5/a to the scale [51]. 
Since we can constrain Sja <0.1 via the size of the magnetic 
moment, then one expects, all things being equal, the trans¬ 
verse scale J± to be suppressed by two orders of magnitude 
relative to the super-exchange contributions to the Ising cou¬ 
pling Since we know 7 ~ 5 K in Dy 2 Ti 207 [13, 16] and 
7 ~ 2 K in Ho 2 Ti 207 [15], this implies 7± < 50 mK. This is 
the main conclusion of our work. 

Estimate from model calculation: The simple scaling ar¬ 
gument above does not include combinatoric factors or other 
unforeseen quirks or cancellations that could favor or disfavor 
the generation of transverse 7± couplings. To check that such 
anomalies do not occur, we carry out an explicit calculation 
to verify this scaling estimate. We work within the charging 
approximation as encapsulated in Eq. (12), considering only 
a single super-exchange path. 

The shortest path between the rare-earth ions and oxygen 
passes through the O' site situated in the center of each tetra¬ 
hedron. Within the Slater-Koster approximation [52], the hop¬ 
ping matrices fi, t 2 can be expressed as f, = r] to where /?, is a 
rotation of the / and p orbitals that takes a set of axes aligned 
to the i local axes into the global frame. The matrices to de¬ 
fine overlaps in the frame aligned along the bond axis and take 
the simple form [toUm’ = 5mm'(<5|m|=ifp/7( + dm=otpfo-)- We note 
that tpf„ is smaller in magnitude than tpfa- and of opposite sign 
[53], though their precise values will not be important to our 
discussion. Since to is diagonal and Rj is unitary, the products 
that appear in Eq. (12) are 






(14) 


as in Ref. [34] . Consideration of further super-exchange paths 
between the other oxygen atoms would change the form of f,, 
but not the overall structure of the interactions. 

To evaluate the super-exchange interactions, we proceed as 
follows: we first construct the CEF ground doublet states of 
Eq. ( 1 ) from the full manifold. Next, we project the one- 
electron operators into this basis. Finally we sum these 
one-electron operators with the interaction constants of Eq. 


(12) evaluated using the f, given above for single bond of the 
lattice. The remaining bonds can be recovered using the lattice 
symmetry. This gives a model of the form shown in Eq. (5). 
At leading order in the d„ coefficients, one has 

Jxx = + Oiad^), 0{a6^), 

7^ = -7 + 0{a5^), -/(/!/ = 0 + 0{a6^), (15) 


where all of the have been considered at the same order in 
the expansion [54]. The super-exchange energy scale 7 is 


7 = 


pfo- 


(U} + A)2 


1 




u: 


2Uj + Up+ 2A 


2187 


121 +96(—) +45o(—) 

xhfo') xhfo-/ , 


(16) 


From Eqs. (15) and (16), we see that our naive scaling argu¬ 
ment of 7± ~ (dlafj does in fact hold in this more detailed 
calculation. We also note that the sign of the super-exchange 
contribution to is negative, as required for Dy 2 Ti 207 
[13, 16]. 

Only the super-exchange contribution to 7^^ has been com¬ 
puted here; when rotating to eliminate the J^z coupling, we 
must also include the large contribution D to coming from 
the nearest-neighbor part of the MDD interactions. We thus 
shift Jzz ^ Jzz + D where D is given in Eq. (8). Using the 
approximate rotations of Eq. (7), we find Jy ~ 0 and 


Jx 


4DJ 


— D — J + 


47^ 

D-J\a) 


(17) 


From the constraint on the moment size, we find 6^la^ < 
0.005 so then 7± < 56 mK. At the temperature scale of in¬ 
terest, monopole excitations [10] are suppressed and quantum 
effects can proceed only through quantum tunnelling within 
the spin ice manifold [55]. The strength of tunnelling g ap¬ 
pears perturbatively at third-order in the transverse coupling 
as gf ~ 1272/7 j [43, 55, 56]. Using the above estimate, these 
quantum effects would only become relevant for temperatures 
< 0.14 mK in Dy 2 Ti 207 . In fact, this estimate could be re¬ 
duced significantly depending on how close the sub-leading 
spectral components of |±) are to their maximal values al¬ 
lowed by the moment size. For example, using the CEF pa¬ 
rameters of Bertin et al. [42], one has 6s ~ 10“^ which is 
much smaller than the bound implied by < 0.025/n from 
the moment, giving a significantly smaller transverse coupling 
7±. Following the same methods, the estimate for tunnelling 
in Ho 2 Ti 207 is also <0.1 mK [41]. 

Outlook: We have presented a general argument as to why 
transverse exchanges in the canonical spin ices are small. 
In particular, we have shown that the generation of the re¬ 
quired high rank multipolar couplings needed to link the 
nearly maximal |+J) states of the CEF doublets are strongly 
suppressed. Using an approximate treatment of the oxygen 
mediated super-exchange interaction, we have provided some 












5 


heuristic bounds on this size of these couplings. The pres¬ 
ence of all these factors can provide an explanation of spin 
ice behavior in the related germanates Ho 2 Ge 207 [32, 39] and 
Dy 2 Ge 207 [29, 39] or the stannates Ho 2 Sn 207 [28, 37] and 
Dy 2 Sn 207 [28, 38]. These criteria also suggest what features 
to look for to move away from the classical spin ice limit [5]. 
For example, compounds with Kramers ions and large Ising- 
like moments, but significant sub-leading components in their 
ground state CEF doublet could present quantum behavior at 
more experimentally relevant temperature scales [5]. In this 
context, the spinel CdEr 2 Se 4 [57] with Ising Er^^ moment is 
an intriguing example where this logic may apply. 
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TRANSVERSE INTERACTIONS IN HOLMIUM TITANATE 


Ground state doublet 

In the case of the crystal field (CEF) ground state is a non-Kramers doublet built from 
the manifold of the 4/^° configuration [1, 2]. The first excited CEF level is ~ 20 meV above 
the ground doublet [2] . This non-Kramers doublet transforms in the Eg representation of the D^d 
site symmetry group and can be written 

1+) = cr|+8) + di|+5) + d2l±2) + 53l+l) + d4l+4) + d5l+V) (1) 

where \M) = |8,M) are eigenstates of the and E operators of the manifold. This doublet is 
naturally Ising-like with 

<cr| E|cr') = 0, <cr| E|cr') = AcrSo-o-', (2) 

where cr = + and the d-factor is 

A = 8-3(6^1+26^ + 36l + A6l + 5d^). (3) 

Including the Ho^'^ Eande factor, ^ = 5/4, the magnetic moment // = is found to be very 
close to the maximal 10 /zb experimentally [1, 3, 4]. This constrains the spectral content of the |+) 
states in Eq. (1) to be predominantly |+8). We can also appeal to more detailed information on the 
CEF structure. Fits to the CEF levels obtained via neutron scattering have yielded estimates of the 
sub-leading components For example, work by Berlin et al. [5] gives 

di= 0.189, fi2 = 0.014, d3 = 0.070, d4 = 0.031, dj = 0.005. (4) 

Rather than rely on these parameters, we shall consider the bound implied by a deviation of the 
moment from maximal by some amount e. In Dy 2 Ti 207 , we used e = O.OI due to the +0.I//B error 
bar from Ref. [6]. For Ho 2 Ti 207 this relation is 6^ < 8el{3n), but given the lack of stated bound 
on the moment, we use a more conservative e = 0.05 which implies 6^ < 0.133/n. 

Effective pseudo-spin 1/2 model 

We proceed similarly to the Dy 2 Ti 207 case in the main text. Since the CEF energy scale is large, 
it is sufficient to first consider only a bare projection of the microscopic Ho^'^-Ho^'^ interactions 
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acting within the full ^Ig manifold into the ground doublets. Sueh a model ean be expressed in 
terms of the pseudo-spin 1/2 operators 

5* = |+)<+|, S^ = ^(|+)<+|-|-)<-|). (5) 

aeting within the low-energy Hilbert spaee spanned by the |+) doublets. The symmetry of py- 
roehlore lattiee eonstrains the form of any sueh projeeted Hamiltonian. For a non-Kramers dou¬ 
blet, a nearest-neighbor model ean be shown to be restrieted to the form [7, 8 ] 

^ [45/55 - 7, (5/57 + 5-57) + 7,, (70-5757+r;;^r^ 7 )l • 

(ij) 

where the bond dependent phases 7,7 are defined in Ref. [ 8 ] . Onee dipolar interaetion are ineluded, 
we will find that 7^^ is the dominant bare eoupling and that it is positive. We will eharaeterize 
deviations from the Ising limit using the nearest-neighbor transverse scale 7+, as quantum effeets 
from 7++ appear at higher order in peturbation theory [ 8 , 9]. Our goal is to estimate the size of 7+ 
whieh is responsible for the quantum dynamies within the low-energy pseudo-spin 1/2 model. 


Ising interactions 


As in Dy 2 Ti 207 , the long-range dipole-dipole interaetions play an important role, In Ho 2 Ti 207 
these take the same form 



(nj ■ zi) (hj ■ z;)], 


(7) 


where r„„ is the nearest-neighbor distanee between the rare-earth atoms, = rt - rj is the dis¬ 
placement veetor between the two sites and z, is the quantization axis along the local cubic [ 111 ] 
direction at site The strength of the eoupling is eharaeterized by ^ The 

largest pieee is the nearest-neighbor eontribution, with coeffieient D = 5^/3 which is ~ 9.4 K in 
Ho 2 Ti 207 [10], slightly larger than in Dy 2 Ti 207 . As in Dy 2 Ti 207 , to obtain a quantitative agree¬ 
ment with experiments, eorreetions to the Ising eouplings must be ineluded. The nearest-neighbor 
eorreetion 7 is the largest, estimated to be ~ 2.08 K [11] for Ho 2 Ti 207 . 


Model calculation 

We follow the same proeedure as in the main text for Dy 2 Ti 207 , this time working with the 
states defined in Eq. (1) drawn from the manifold appropriate for Ho^'^. Aside from the 
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construction of these states, the methods are identieal. We find 



f^) +Oia6\ 

(8a) 

\ 1 

\ a / 


1 + 

1 + 

II 

1 

1 + 

1 + 

l(^) +Oiad\ 

(8b) 

\ hfo- / 

'\a / 


7^2 = —J + O(cr^d), 


(8c) 


where all of the have been eonsidered at the same order in the expansion. The super-exehange 
energy seale J is given by 


J = 


pf<^ 

(f/; + 


1 




+ 




1 


2187 


121+48 


/ 

\ hfo- 1 


+ 207 


\4\ 


I 

^Pfo-j J 


(9) 


The p+ and p++ are rational funetions of the ratio v = tpfj^jtpfa- defined as 


P±(+) = 


9x^(160 - 8x + Ix^) 
8(121 +48x2 + 207x4)’ 


P±±(+) = 


9x^(68 - 16x - 7x^) 
16(121 +48x2 + 207x4)' 


( 10 ) 


Reasonable values of the Slater-Koster ratio satisfy -1 < tpfnitpfa- < 0. For this range, we ean 
bound these funetions as 


0 < P±{tpf„ltpfa-) < 0.6, 0 < P±±{.tpfnltpfa-) < 0.15. (11) 

As in Dy 2 Ti 207 , we see that sealing from the main text, J± ~ ( 6 /a)^J holds with no large eombi- 
natorie faetors or spurious caneellations. We also note that the sign of the super-exehange contri¬ 
bution to the is again, eorreetly negative [10]. Setting the super-exehange seale to 7 = 2.08 K 
from the experimental eonstraints [10], we then ean estimate 7+ and 7++ as 

7± < 1.25 , |7±J<0.31K-(^) . (12) 

For a moment size that deviates from maximal by 6 = 0.05, this ratio is bounded as { 65 /a)^ < 0.027 
giving 

7± < 34 mK, |7±±| < 8.5 mK. (13) 

Here the relevant transverse seale is the 7+ that appears in the nearest-neighbor model (6). With 
7^ = D - 7 ~ 7.32 K, the implied tunnelling seale ^ is 0.1 mK and thus is not experimentally 
relevant. We note that these estimates for 7+ and 7++ seale linearly in the deviation from the 
maximal moment. For the CEF parameters of Ref. [5] we have {65/ a)^ ~ 10“^, and one would 
have to look at next order in the 6 „ to get at the true leading behavior. It is unelear if the fitting of 
the CEF energy levels is suffieiently preeise to provide a useful estimate of the true size of 65. 
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GENERAL DISCUSSION OF SUPER-EXCHANGE 


We use the notation of the main text. As diseussed there, we eonsider the bare Hamiltonian of 
the two rare-earth sites and the intermediate oxygen atom 




(14) 


We perturb this with the hybridization terms 


^ = Z + KapflPp + tiapPlfip + h,apflPp\ , (15) 

a/} 

Exehange interaetions between the rare-earth ions eome in at fourth-order in V. We will follow 
Lindgren [12] and write the effeetive Hamiltonian using the ground state projeetor P and resolvent 
operator R 


P = ^\a) {a\, 

aeO 



(16) 


where lor) are ground states of Hq with energy Eq and \a) are the exeited states of Hq with energies 
Ea- The fourth-order effeetive Hamiltonian is then [12] 


//eff = PVRVRVRVP - PVR^VPVRVP - PVR^VRVPVP 

- PVRVR^VPVP + PVR^VPVPVP (17) 


Now sinee PVP = 0, we ean drop all but the two terms PVRVRVRVP and PVR^VPVRVP. The 
seeond of these faetors as {PVR^VP){PVRVP) and serves to eaneel non-extensive eontributions 
that arise from the first term. 

We thus ean simply use //eff ~ PVRVRVRVP, taking to eare to eonsider only this set of (ex¬ 
tensive) perturbative proeesses. At lowest order, we need only eonsider the exeited eonfigurations 
that ean be built by moving at most two eleetrons between the rare-earth and oxygen sites. These 
types of intermediate states ean be understood by eonsidering the possible values for the oxygen 
eleetron number. The range of aeeessible states are illustrated in Figures 1, 2 and 3. Starting 
from the ground state eonfiguration, one follows the eolored arrows to traee out a given 

exehange path. Spin and orbital indiees on the / and p eleetron operators have been suppressed for 
elarity. We see only two distinet types of proeesses are needed at fourth order, proeess 1 (Figure 
1 ) whieh ereates at most one hole on the oxygen and proeesses 2 and 3 (Figure 2 and Figure 3) 
whieh ereate an intermediate state with two holes on the oxygen. A seeond set of equivalent paths 
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FIG. 1: Process 1 involves intermediate states with only a single hole on the oxygen site. Scales 
roughly as ~ PiUf + Apf)~^(2Uf)~^ 


ean be obtained by switehing the role of sites 1 and 2 in the proeesses. Process 1 corresponds to 
the pieees of //eff of the form 

P [pU\f,] R [f^t,p] R [p^lf,] R [fjt,p] P. (18) 


Each appearance of the resolvent operator R ean be simplified given the restrictions created by 
the hybridization operator V. For example, the piece P can only involve states of the form 

sinee the ground state manifold contains only states of the form /"p^/”. When aeting on 
these states, the resolvent is simplified to an operator that acts only on the states of the rare-earth 
at site 1 


RfA = Yj 

a6/"+‘ 


\a)i {a\i 

Ef,a - E/,0 + Ap/ 


(19) 


where \a) are states of the eonfiguration and £/,« - E/,o are their energies relative to the /" 
ground state, and Apf is the atomie potential of the oxygen site. An identieal operator P /2 can be 
defined that acts only on site 2. The innermost resolvent can be simplified in the same way. Here, 
both rare-earth sites are involved, but the oxygen state is not. It can be replaeed by the operator 


E E 


|q)i (q|i \ a')2 ( a'\2 
Ef,a + Ef\a' - '^Ef^o 


( 20 ) 


This allows us to write the process 1 contribution as 


P 



( 21 ) 


Now since the resolvents no longer depend on the oxygen degrees of freedom we ean regroup this 
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FIG. 2: Process 2 involves intermediate states with two holes on the oxygen site. Seales roughly 
as ~ t\Uf + Apfr\2Uf + 2Apf + UpT^ 

into oxygen and rare-earth parts 

E E [p.pipfpypp] iPfA^Rfjfi>^R‘;}fvPnfiPf] I'U, (22) 

appv a'j3'p'v' 

where we have broken up the projeetor P into rare-earth and oxygen parts as P = PpPf- The 
oxygen part is trivial with Ppplp/jflpvPp = Spydap- The final eontribution from proeess 1 is then 
given by 


Pf 


E /i [vT ['I'T" 


yafSpv 


/• 


(23) 


We earry out the same argument for proeess 2. Following the arrows in Figure 2 we find the 
eontribution is 

P [p^lf,] R [pU\f,] R [f^t,p] R [f,t,p] P (24) 

The outermost resolvents ean simplified in the way. The innermost resolvent needs a new definition 

\a)i (all |a')2 (a'|2 


n(2) _ 
^// = 


Z Z 


Qg/n+l Q'gyi 


Ef^a + Ef'a' - 2P/,o + 2Apf + Up 


(25) 


to take into aeeount the additional oxygen holes and Coulomb repulsion in the y”+ip 4 yn+i §^^^£ 3 ^ 
as diseussed in the main text. This then gives the operator 


p[p^Afi\Pf 2 [p^^lfi]pff [flhp]Pf,y [flhp] P- 

Onee again we ean faetor out the oxygen part as 


(26) 


'y'l ^ l^PpPaP^sPfjPvPpl [Pffla'PplfljS'P^ffflp'Pp^fL'Pf^ [ 4 ,aQ'Tl,/ 3 ygT 2 ,i/'^^l,v'v] • ( 22 ) 
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FIG. 3: Process 3 involves intermediate states with only a single hole on the oxygen site. Seales 
roughly as ~ t\Uf + Apf)-H2Uf + 2Apf + Up)-^ 


The oxygen part is a slightly less trivial this time with Ppplp^^PpPyPp = S^pdav - SapS^y We thus 
have 


Pf 


E[‘>41” - ['24 

\aPjiv 


y/J 


Pf 


(28) 


Finally, for process 3, we ean reuse all of the building bloeks that we have put together so far. 
From Figure 3, we directly write the eontribution 


P[p^Afi][p^Af^Pff [fihp]Pf,y [fi hp] P 


(29) 


Following an identieal argument to the previous seetion we arrive at 


Pf 


E /,„«/../ 2 ««iX%' 44 ([' 24 ] ['■'il” - [vlT ['i4] 

\a/}fiv 


y/3 


(30) 


by simply swapping the site 1 and site 2 labels in part of the exehange path. 


MULTIPOLE RANKS AND RESOLVENT OPERATORS 

From the expressions above, we see that the interaetions between the rare-earth ions ean be rep¬ 
resented by operators of the form ~ mediated by the oxygens, intermingled with resolvents 
sueh as Rf i, and As discussed in the main text, we eonsider the case where the resolvents 
do not eontain the CEF potential and are thus alforded the full spherieal symmetry of the free-ion. 
Beeause of this symmetry, we ean write the resolvents in terms of spherieally symmetrie projeetors 
into each set of degenerate levels. Expanding out these sums and re-arranging terms, we find that 
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contributions to each process take the form 


'(/.lO,/,;,) ( 4 , 02 / 2 ,,), (31) 

where 0 \ and O 2 are operators acting on rare-earth sites 1 and 2 built from rank -0 spherically 
symmetric projectors. The same rank eounting outlined in the main text within the eharging ap¬ 
proximation thus holds true; operators of the form involve multipoles of rank-7 or lower. 


CHARGING ENERGY APPROXIMATION 


Here we give the details that link the general results to the final expressions for the eharging 
approximation given in the main text. Reeall that in this approximation we postulate the energies 
of the /" states predominantly depend on the total electron number n. With this approximation we 
write 


Rf.i 


z 


\a)i {a\i 


l-P 
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—. t/7 -I- Apf 

aef"+' f PJ 




(32) 


sinee Ef^a - Ef^ ~ Uj for all states \a) e Sinee the strueture of the perturbation forees us 
to be in the state, we ean extend the sum to all excited states and thus arrive at the projeetor 
1 - Pf. Applying the same approximation to the innermost resolvents, we find 

l-Pf 


R 


( 1 ) ^ / 

// u}^u-/ 


R 


^-Pf 


m __ 

f'f 2Uj + Up + 2Ap/ 


(33) 


where Uj is the energy differenee Ef^a - Ef^ when \a) is a state in the manifold. We ean 
now use these resolvents to simplify the eontributions from the three proeesses (Figs. 1, 2, and 
3) eonsidered earlier, and shown in Eqs. 23, 28, and 30. In eaeh expression the ground state 
projeetors Pf drop out, and we ean write the three parts 


1 : 


1 


Pf 


\a/jfjv 
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3 : 


(U} + ApfnUj + Uj) 

1 

(U} + Apfn2U} + Up + 2Apf) 
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Pf 
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Pf 


We ean further massage the expressions for eaeh proeess to put the / operators into a standard 
order. Finally, the swapping of 1 and 2 generates an identieal contribution and thus a factor of 
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two. Collecting the contributions from all three processes and dropping any single-site operators 
generated by commuting the pieces, we have the final result for //eff 


aPfiv 


(34) 


where the interaction matrix is defined as 
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■apiiv 
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([/; + 


u}^u- 
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2Ut + U„ + 2A 
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•pf 


] [. [. d]“" ^ 

hA 


t,t\ 

ap - 

jpiMj ri^2j 2Uj + Up 

+ 2Apf 


(35) 


as given in the main text. An alternative route to the super-exchange interactions would be to in¬ 
tegrate out the oxygen degrees of freedom first. This down-folding gives an effective /-/ hopping 
of order = t^lApf. Within the approach discussed here, this corresponds to taking Apf » Uj, 
eliminating all but process 1 at leading order. This then takes the form of a second-order result in 
the effective hopping t^fi 


OTHER MICROSCOPIC INTERACTIONS 

Here we outline some other possible sources of interactions with the J = 15/2 and J = 8 
manifolds of Dy^”^ and Ho^"^. 

Magnetic multipole interactions 

Being of magnetostatic origin, the scale of the magnetic multipole interactions (MMI) are sim¬ 
ple to address. These are only directly relevant for Dy 2 Ti 207 , as the transverse couplings in 
non-Kramers doublet relevant for Ho 2 Ti 207 are time-reversal even. Generally, the scale of the 
interaction between a rank-A^ and rank-^' magnetic multipole is approximately given by [13] 

Mkk' (^ 6 ) 

where is the order K Stevens’ parameter [14], (r^) is an atomic radial integral [15], ju. is the 
magnetic moment of the ion and R is the distance between the ions. This expression is somewhat 
conservative as we have assumed the multipole moments are maximal, including factors of and 

in where J = 15/2. Projecting into the ground state doublet will quench these moments 
and further reduce the estimate. This is especially true for Dy 2 Ti 207 where the ground state 


10 











doublet is primarily composed of 1+15/2) states and thus has predominantly rank-15 moments 
that could contribute the transverse couplings. Using ~ 10 jUb, and we find these MMI are small 
beyond the dipole-dipole term Mu ~ 1.4 K. The more useful number is is D = 4 (5 A1ii /3) ~ 9 K 
as discussed in the main text. For dipole-octupole coupling one finds AI 13 ~ -16 mK while 
for octupole-octupole one finds AI 33 ~ 0.17 mK. Given the factor of higher-rank 

interactions are quickly suppressed. The nearest-neighbor dipole-octupole coupling is largest and 
necessarily projects into exchanges of the form ~ 5 y + 5 Finally, one must take into account 
that when the octupole moments are projected into the ground state doublet, their contribution is 
further reduced. The suppression factors are ~ d/cr, which reduce the dipole-octupole terms by 
a factor of ~ 0.1 and the octupole-octupole terms by ~ 0.01. The induced transverse terms after 
rotation is thus of order < 10“^ mK. The quantum effects are then completely negligible, and thus 
the MMI can be safely ignored beyond dipole-dipole interactions discussed in the main text. 

Electric multipole interactions 

Electric multipole interactions (EMI) are only directly relevant for Ho 2 Ti 207 as the transverse 
couplings in the Kramers doublet relevant for Dy 2 Ti 207 are time-reversal odd. Generally, the scale 
of the EMI between a rank-K and rank-K' electric multipole is approximately given by [16] 

&KK' ~ XKXK'{r^){r^ ) [^1 ’ ( 37 ) 

where e is the elementary charge, 60 is permittivity of vacuum and the remainder of the notation is 
the same as in the previous section. We have included factors of and to account for possible 
matrix element effects as conservatively as possible. The largest of these couplings is the electric 
quadrupole-quadrupole coupling at order &21 ~ 63 mK. Due to the composition of the ground 
doublet, transverse couplings originating from EMI will carry additional factors of ~ (d/or)^, as 
in the MMI. The scale of these contributions are then < 1 mK. The higher order interactions are 
even smaller and so the electric multipolar interactions are not a significant source of transverse 
couplings in Ho 2 Ti 207 . 

Lattice mediated interactions 

Additional interactions could be mediated through spin-phonon coupling [17]. Due to time- 
reversal invariance, these can only directly generate interactions between multipoles of even rank. 
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For Dy 2 Ti 207 these vanish under projection into the ground state CEF doublet, and thus do not 
contribute directly to the transverse terms. For Ho 2 Ti 207 , such even-rank interaction can con¬ 
tribute directly to the transverse terms. To evaluate the plausibility of such a scenario, we consider 
two models of spin-phonon interaction, a site-phonon and bond-phonon model. In the site-phonon 
model, we consider only the coupling of the on-site multipoles to the local lattice distortion. In the 
same way as the CEF, when represented as Stevens’ operator equivalents, these can only include 
operators up to rank-6. When the site-phonon is integrated out, these operators are squared so only 
interactions of rank-12 or less can be generated. In a bond-phonon type model, we imagine the 
lattice distortion modifying the exchanges already present in the system. Given the argument in 
the main text, we expect these to be predominantly rank-7 or lower and thus when the phonon is 
integrated out the induced interactions are rank-14 or smaller. Since these are not the direct rank- 
16 interactions needed to connect |+8) to |+8), they are suppressed by the sub-leading doublet 
components by same mechanisms discussed in the main text. Additionally, there will be a further 
suppression from the energy scale of the phonon that was integrated out. We thus can safely ignore 
the spin-lattice couplings in our analysis. 


Direct exchange 

We have accounted for some of the inter-site Coulomb interaction in our treatment of EMI. 
Other contributions such as direct exchange between the rare-earth ions is likely to be of the same 
order or smaller given their large separation. Generally, one could write pair-wise direct exchange 
contributions in the form 


E E ( 38 ) 

ij aPfiv 

As we have seen, the one-electron form of the operators at each site restrict this to generate multi¬ 
polar interactions of rank-7 or less. Thus we have the same suppression factor (d/cr)^ as discussed 
in the main text. Given the scale of the EMI, we then expect these terms to be negligible. We 
note that direct exchange interactions between the rare-earth and oxygen sites could be present. 
However, due to the need to include additional electron transfers to generate interactions and the 
filled shell of the oxygen site, we would not expect these terms to be important. 
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